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By using Supersymmetric Quantum Mechanics and Semiclassical Quantization, one may argue 
that the low-lying excited states of any quantum system can be modeled by a set of harmonic 
oscillators. In the present paper, we fit the experimental excitation spectra of atoms with atomic 
number 2 < Z < 36 to a simple harmonic oscillator model with two parameters: the number of 
degrees of freedom, d, and the effective frequency, oj. The obtained hu takes values around 0.03 
(in atomic units), whereas d shows clear shell filling effects, that is, takes high values for the noble 
gases, suggesting collective oscillations of the electrons occupying the last shell. 



I. INTRODUCTION 

El presente trabajo tiene la intention de ser un mo- 
dest o homenaje al Prof. Marcos Moshinsky, fallecido re- 
cientemente, quien fue un defensor del papel jugado por 
el oscilador armonico en la Fisica moderna y un pionero 
en la construction de bases de funciones de tipo oscila- 
dor para la description mecanocuantica de los nucleos 
atomicos [1 . 

En contraste con el rigor matematico del Prof. Mos- 
hinsky, vamos a seguir un estilo muy cualitativo e inten- 
tar, con argumentos simples y ajustes a datos experimen- 
tales, dar respuesta a la interrogante de si el espectro de 
excitaciones de los atomos se puede modelar con oscila- 
dores armonicos. Est a interrogante se puede desdoblar, 
en realidad, en dos aspectos muy interesantes, que abor- 
daremos someramente. 

Por un lado, es bastante extendida en la comunidad 
de ffsicos la intuition de que cualquiera sea el sistema 
cuantico o el potencial de interaccion entre particulas, 
el espectro de excitaciones a bajas energias es siempre 
el de un sistema de osciladores. Daremos un argumento 
a favor de este cuadro cualitativo apelando a la cuanti- 
zacion semiclasica en un potencial efectivo sugerido por 
la denominada Mecanica Cuantica Supersimetrica. Es- 
te argumento fue formulado en [2 y aplicado a sistemas 
modelos de tres particulas. 




Figura 1: El Prof. Marcos Moshinsky, premio Principe de As- 
turias de 1988 y Premio UNESCO de las Ciencias en 1997. 



Por otro lado, es tambien muy atractiva la idea de 
ajustar con pocos parametros efectivos datos experimen- 
tales o calculados de una magnitud fisica. Esto es muy 
comun en el quehacer de la Fisica experimental. Como 
ejemplo en el ambito de la teoria de moleculas, muy cer- 
cano al espiritu del presente trabajo, podemos citar el 
hecho conocido de que el espectro de oscilaciones de las 
protemas puede ajustarse muy bien con una sola cons- 
tante de fuerza para describir la interaccion entre cual- 
quier par de atomos [3]. En nuestro trabajo utilizaremos 
la excelente compilation del NIST sobre los espectros de 
excitaciones de atomos [4 y ajustaremos cada espectro 
con solo dos parametros: la frecuencia (uj) y la dimension 
efectiva (d) del oscilador. Para un atomo como el Na, por 
ejemplo, con un solo electron de Valencia, uno esperaria 
d ~ 3 y co describina la interaccion de ese electron con el 
resto de los electrones internos. En el caso mas general, 
d/3 seria un indicador del mimero efectivo de electrones 
que participan en las oscilaciones. 



A. El cuadro de osciladores para los estados 
excitados 

Consideremos un sistema cuantico de varias particulas, 
descrito por el Hamiltoniano: 

donde V da cuenta de la interaccion. Si la funcion de 
onda del est ado base del sistema, ipo = exp—W/h y la 
energia, Eq, son conocidas, podemos escribir la siguiente 
expresion para el potencial: 

v ~ E » = E or { ( w)2 - h ^ w )) - ( 2 ) 

Esta expresion se conoce como "representacion del esta- 
do base" y tiene su origen en la denominada "Mecanica 
Cuantica Supersimetrica" [5], la cual se construye por 
analogia con la super simetria en la Fisica de particulas, 
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pero en realidad trata sobre relaciones de casi isoespec- 
tralidad entre hamiltonianos. El denominado "superpo- 
tencial" : 

£/ = E^( v ^) 2 > (3) 

i 

tiene las siguientes propiedades: 

- U > aiin cuando V no sea acotado por debajo. 

- El mmimo, U = 0, corresponde a ViW = 0, es de- 
cir la configuracion de maxima probabilidad. Esto 
es una extension del concepto de radio de Bohr a 
sistemas con muchos grados de libertad. 

- U difiere de V — Eo en terminos que son formalmen- 
te de orden superior en H. Lo que sugiere que la 
regla de cuantizacion de Einstein-Keller-Maslov: 

I y^Pi • dn = (n a + |m a )fi (4) 

podrfa aplicarse en U. 

- El est ado base, al que corresponde energfa de exci- 
tation igual a cero, es descrito exactamente por la 
Ec. Q con indices n a = 0, m a = 0, por eso dicha 
condition de cuantizacion es buena no solo cuando 
n a es grande (como sucede con toda semiclasica) , 
sino tambien cuando los n a son pequehos. 

- La Ec. Q con m a = es exacta en toda una clase 
de problemas de una particula en una dimension, 
que incluye al potencial de Coulomb. 



n 2 




Figura 2: Estados excitados de un oscilador bidimensional. 
Los estados con energfa menor o igual que Ai? conforman un 
triangulo, cuya area viene dada por la parte derecha de la Ec. 

En primera aproximacion, las frecuencias de pequehas 
oscilaciones alrededor del mmimo de U dan el espectro 



de excitaciones del sistema. En vez de seguir el procedi- 
miento de calcular la funcion de onda, construir U, etc, 
nosotros tomaremos el espectro medido del atomo y lo 
ajustaremos con dos parametros: una frecuencia prome- 
dio, ou y una dimension efectiva, d, que da cuenta del 
mimero de grados de libertad que participan de las osci- 
laciones. 



B. El oscilador en d dimensiones 

Consideremos un oscilador de frecuencia uj en d dimen- 
siones. La constante h se toma igual a uno. El mimero de 
estados con energfa de excitacion menor o igual a AE es: 

^ (AE/u,) d 

iy estados ~ ^ • \ D ) 

En efecto, AE/uj es el mimero maximo de cuantos de los 
estados inclufdos y (AE/uS) d jd\ el volumen de la region 
con energfa menor o igual a AE en el espacio de los mime- 
ros cuanticos n a . En d = 2, por ejemplo, la parte derecha 
de la Ec. ([5| da el area del triangulo correspondiente en 
la Fig. [2) 




o!o o!l 0,2 t 3 OA o!s 0,6 
Energfa (u.a.) 



Figura 3: Espectro de excitaciones del Ar. Arriba: ajuste ini- 
cial del parametro d con las excitaciones mas energeticas. De- 
bajo: ajuste simultaneo de los parametros d y uj para todo el 
espectro. 

Al ajustar un espectro atomico permitiremos al 
parametro d tomar valores no enteros, lo cual puede en- 
tenderse como un valor promedio entre las distintas d que 
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aparecen a distintas energias de excitation. Por eso en el 
denominador de la Ec. ([5| sustituiremos d\ por T(d + 1). 

II. LOS ESPECTROS DE EXCITACIONES DE 
ATOMOS 

Se utilizaron los espectros de niveles de energia ofreci- 
dos por el NIST [4] para los atomos de la Tabla Periodica 
con 2 < Z < 36, desde el He hasta el Kr. Para cada ato- 
mo tomamos los primeros 30 niveles con el objetivo de 
ajustar el espectro segun la formula En el numero 
de estados con energia menor o igual a AE se incluye 
explicit anient e la multiplicidad de cada estado debida al 
moment o total J. 

El numero 30 de niveles incluidos es completamente 
arbitrario y, en atomos distintos, puede corresponder a 
diferentes proporciones de excitaciones intracapas e in- 
tercapas, las cuales exhiben diferentes comport amientos 
en la energia de excitation [6]. Sin embargo, la experien- 
cia nos sugiere [7 fijar el numero de niveles al comparar 
sistemas cuanticos con distintos numeros de partfculas. 

En la Fig. [3]se muestra, en el ejemplo del Ar, el proce- 
dimiento seguido para determinar los parametros d y uo. 
Recordar que nos limitamos a los primeros 30 niveles de 
energia. Utilizando inicialmente escala doble logaritmica 
y fijandonos solo en las excitaciones de mayor energia ha- 
llamos un valor para d. Ese valor initial es usado como 
punto de partida para el ajuste de todo el espectro con 
la formula (|5|. Los parametros obtenidos se muestran en 
la Tabla H 
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Figura 4: Espectros de excitaciones de bajas energias y ajustes 
segun la Ec. j5} para diferentes atomos. 

En la Fig. [4] most ramos que la calidad del ajuste varia 
para diferentes atomos. En atomos con capas abiertas, 
las excitaciones intracapas no son muy bien descritas por 
la ley Si nos hubieramos extendido hasta mas al- 
tas energias tomando, por ejemplo, 100 niveles en vez de 
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Z 


d 


w 


Pot. Ioniz. ( 


He 


2 


28.19 


0.066 


0.9040 


Li 


3 


8.21 


0.027 


0.1982 


Be 


4 


7.33 


0.048 


0.3428 


B 


5 


6.25 


0.044 


0.3051 


C 


6 


3.94 


0.048 


0.4140 


N 


7 


3.96 


0.062 


0.5344 


O 


8 


10.32 


0.062 


0.5007 


F 


9 


9.36 


0.082 


0.6406 


Ne 


10 


13.81 


0.089 


0.7929 


Na 


11 


9.81 


0.024 


0.1890 


Mg 


12 


8.02 


0.038 


0.2811 


Al 


13 


6.90 


0.030 


0.2201 


Si 


14 


8.38 


0.034 


0.2997 


P 


15 


8.90 


0.045 


0.3856 


s 


16 


10.15 


0.045 


0.3809 


CI 


17 


11.61 


0.053 


0.4768 


Ar 


18 


18.56 


0.052 


0.5794 


K 


19 


8.46 


0.021 


0.1596 


Ca 


20 


5.71 


0.027 


0.2248 


Sc 


21 


3.43 


0.010 


0.2413 


Ti 


22 


1.73 


0.003 


0.2511 


V 


23 


5.46 


0.010 


0.2481 


Cr 


24 


6.57 


0.016 


0.2488 


Mn 


25 


2.88 


0.013 


0.2733 


Fe 


26 


6.08 


0.013 


0.2906 


Co 


27 


1.61 


0.004 


0.2898 


Ni 


28 


7.60 


0.019 


0.2809 


Cu 


29 


3.64 


0.030 


0.2841 


Zn 


30 


10.74 


0.043 


0.3454 


Ga 


31 


8.08 


0.030 


0.2206 


Ge 


32 


9.37 


0.033 


0.2904 


As 


33 


8.68 


0.042 


0.3599 


Se 


34 


7.13 


0.046 


0.3586 


Br 


35 


14.11 


0.043 


0.4344 


Kr 


36 


12.25 


0.059 


0.5148 



u.a.) 



Cuadro I: Parametros ajustados d y uo y potenciales de ioni- 
zation medidos para los atomos con 2 < Z < 32. 



30, los ajustes hubieran parecido cualitativamente mucho 
mejores. 

La Fig.|5]ilustra la dependencia del parametro d con el 
numero atomico Z, evidenciando claros efectos de llenado 
de capas. Los atomos con capas llenas tienen valores de d 
mas altos. Si interpretamos d/3 como el numero efectivo 
de electrones que participan en las oscilaciones, vemos 
que para el He este numero es 9 (un valor no ffsico), 
pero para el Ne es casi 5, para el Ar es 6 y para el Kr 
es 4. Todos estos son numeros razonables que indican 
que en las excitaciones de mas bajas energias participan 
colectivamente los electrones de la ultima capa. 

Por otro lado, el parametro cu es la curvatura efecti- 
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Figura 5: Relacion entre la dimension efectiva d y el numero 
atomico Z. 



va del superpotencial alrededor del mmimo U = 0. Es 
mas dificil relacionar este parametro con alguna magni- 
tud atomica. En la Fig. [6] comparamos el potencial de 
ionizacion del atomo con el valor ajustado de uj. Estas 
magnitudes deben estar relacionadas, aunque el poten- 
cial de ionizacion corresponde al tope del espectro discre- 
to, es decir a mas altas energias de excitacion. La curva 
contmua en la Fig. es solo una guia visual, pero muestra 
que existe cierta correlacion entre las dos magnitudes. 




0,02 0,04 0,06 

Frecuencia de Oscilacion 



Figura 6: Relacion entre el Potencial de Ionizacion y la fre- 
cuencia de oscilacion u>. La lmea roja es solo una gufa visual. 



III. CONCLUSIONES 



El presente trabajo se ha nutrido de varios ingredien- 
tes. Por un lado, utilizando la denominada Mecanica 
Cuantica Supersimetrica, argumentamos que las excita- 
ciones de mas bajas energias de cualquier sistema cuanti- 
co corresponden a osciladores armonicos. Y planteamos 
un modelo sencillo para el numero de estados excitados 
con energia de excitacion menor que AE. Este modelo 
tiene dos parametros: el numero de grados de libertad, d 
y la frecuencia promedio de los osciladores, uj. 

Por otro lado, utilizando los datos experiment ales de 
los espectros atomicos ofrecidos por el NIST, determi- 
nants los parametros d y uj para todos los atomos con 
2 < Z < 36. Para ello se utilizaron los primeros 30 ni- 
veles excitados de cada atomo. En realidad, si hubiera- 
mos utilizado mas niveles, por ejemplo 100, la calidad 
del ajuste hubiera sido muy superior. No lo hicimos por- 
que no tenemos respuesta a la pregunta de por que apa- 
rentemente no hay efectos anarmonicos a energias mas 
altas. Es decir, por que los atomos se comportan casi 
como osciladores perfectos. Las variaciones del parame- 
tro d se relacionan con el llenado de las capas atomicas y 
"muestran" el caracter colectivo de las excitaciones en los 
atomos con capas llenas. El parametro por su parte, 
evidencia cierta correlacion con el potencial de ionizacion 
del atomo. 

Los ajustes de los espectros atomicos con el modelo 
simpliflcado de osciladores, dado por la Ec. (|5|, son cuali- 
tativamente buenos. Lo cual trae de vuelta, como muchas 
veces sucede en la Ffsica y en otras ciencias, a un viejo 
modelo, descartado en su momento, el modelo de Thom- 
son para los atomos, donde el potencial de confinamiento 
era armonico. 
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